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Abstract 



A 2n-dimensional Lax integrable system is proposed by a set of specific spec- 
tral problems. It contains Takasaki equations, the self-dual Yang-Mills equations 
and its integrable hierarchy as examples. An explicit formulation of Darboux 
transformations is established for this Lax integrable system. The Vandermonde 
\q ■ and generalized Cauchy determinant formulas lead to a description for deriving 

0^ | explicit solutions and thus some rational and analytic solutions are obtained. 



1 Introduction 



Darboux transformations provides us a purely algebraic, powerful method to construct 
solutions for systems of nonlinear equations Jl] • Its key is to expose a kind of covariant 
properties that the corresponding spectral problems possess. There have been many 
tricks to do this for getting explicit solutions to various soliton equations including 
KdV equation, KP equation, Davey-Stewartson equation, Veselov-Novikov equation 
etc. (see for example 0,0,0,0,0 and references therein). Darboux transformations 
can also be applied to generating multi-soliton solutions to soliton equations and the 
Darboux covariance makes it possible to construct a series of exactly solvable systems 
of supersymmetric quantum mechanics Q. 

In this letter, we would like to establish a kind of Darboux transformations for a Lax 
integrable system in 2n-dimensions which we will introduce. The dimension reductions 
of this system contain some interesting and important equations as examples. Among 
them are Takasaki equations J7|, the self-dual Yang-Mills (SDYM) equations, and the 
self-dual Yang-Mills hierarchy || etc. It is well known that the SDYM equations (even 
the SDYM hierarchy) may be reduced to many integrable equations in 1 + 1 dimensions 
and in 1+2 dimensions when certain symmetry conditions are imposed (see for instance 
0, 0). Therefore our Lax integrable system also includes a lot of integrable soliton 
equations in 1 + 1 dimensions and in 1 + 2 dimensions. Recently there is still an 
considerable interest in the aspect of symmetry reductions of the SDYM equations 



1 



sec 



10 , [11], a quite detailed list of the relevant references is included in fl0"|). This 



also increases, to a great extent, the validity of the Ward's conjecture [12|: many (and 
perhaps all?) integrable equations may be derived from the SDYM gauge field equations 
or its generations by reduction. 

This letter is organized as follows. In Section 2, we derive a Lax integrable system 
starting from a set of specific spectral problems and display a few concrete systems 
of nonlinear equations. The corresponding Darboux transformations are established in 
Section 3 and an explicit description of a broad class of solutions is proposed by means 
of the resulting Darboux transformations. Section 4 contains some further discussions 
and two remarks, along with a more general system whose spectral problems include 
negative powers of a spectral parameter. 



2 Lax integrable system 

Let the differential operators Lj, 1 < i < n, be defended by 

U = Li(X) = _-.,(*)_ = _- (g 0it A'«)^, 1 < i < n, (2.1) 

where the coefficients a ik , < k < M, 1 < i < n, are constants, andp = (p±,P2, • • • ,p n ), 
x = (xi,x 2 , ■ ■ • ,x n ) G(C n or E" are two vectors of independent variables. We consider 
a set of spectral problems 

M 

L,(A)^ = -A(A)# = -(£ Aa\ l )V, 1 < i < n, (2.2) 

1=0 

where A is a spectral parameter, ^ is a N x N matrix of eigenf unctions, and An, 1 < 
i < n, < / < M, are all N x iV matrices of potential functions depending on p, x. 
This spectral problem contains Takasaki case 0: M — 0, aj(A) = A, 1 < i < n; Gu 
and Zhou case [|]: a^X) =0, 1 < i < n, Ai(X) = XJi + Pi, 1 < i < n — 1, where Jj, Pi 
are the matrices given in ||; and Gu case Oj(A) = A, 1 < i < n — 1, a n (A) = 
0, Aj(A) = v4j , 1 < i < n — 1. We assume that n, N > 2 in order to obtain nonlinear 
integrable systems (we shall explain this later). Noting that Li, 1 < % < n, are all 
linear operator, we can calculate that 

M 

L i (A)L i (A)* = J L i (A)(-X)A I i4 il *) 

1=0 

M M 

= -E A%Li(A)* - E A%(A)A0* 

M M M 04 a/, 

= £ A% £ A% fc * - £ A<(^ - *(A)^)* 
/=o fc=o z=o WX J 

= E A ™( E ^,0*-£a^*+£a^( E %^)*- 

m=0 k+l=m m=0 U t ) j m=0 k+l=m udj j 

0<k,KM 0<k,KM 
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Therefore we see that the compatibility conditions 

L j (X)L i (\)^ = Li(X)Lj(X)ty, 1 < i, j < n, (2.3) 

are expressed as 

2M M Q4. 2M a A 

EA m E M^-Ea^+E^ 1 E ^ 

m=0 k+l=m m=0 u Pj m=0 k+l=m ux j 

0<k,l<M 0<k,l<M 

2M M 04 2M 

= Ea™ E V,-EA m ^ + E^ +1 E ^if- (2-4) 

m=0 k+l=m m=0 m=0 k+l=m ux i 

0<k,l<M 0<k,l<M 

Equating coefficients of the terms A ; A m , 1 < m < M; A m , M + 1 < m < 2M 
and \ 2AI+l in the above equation lead to the following Lax integrable system in 2n- 
dimensions: 

r) A r) A 

[40,^1-^ + ^ = 0, l<ij<n, (2.5) 

E [Ak,A Jl }-^ + ^+ e (a iJk ^ - aJ^)= o, 

fc+Z=m ^Pj °W fc+«=m-l 

0<k,l<M 0<k,l<M 

1 < m < M, 1 < i,j < n, (2.6) 
E [A4k,Afl] + E ( a ^^- aifc ^)= ' 

k+l=m k+l=m-l ux j ux i 

0<k,l<M 0<k,l<M 

M + 1 < m < 2M, 1 < i,j < n, (2.7) 

dA iM OA 

ajM ~dx- am dx ' t,J n ' ^ ^ 

where [■, ■] denotes Lie bracket of matrix Lie algebras. Evidently for the case of N = 1, 
the commutators of the matrices A^, Aji are all equal to zero and thus the above 
system is simplified into a linear system, which isn't what we need. Besides, if n = 1, 
the above system holds automatically and thus doesn't need any consideration. The 
Lax integrable system consisting of (|2.5| ), ( p.6|) , ( p.7|) and (|2.8|) includes some interesting 
and important systems of equations as examples. For instance, we have a few systems 
of equations as follows. 

Example 1: M = 0, Gtj(A) = 1, 1 < i < n, corresponds to Takasaki case |7[ which gives 
rise to the equations 

[A i0 ,A j0 ] - + = 0, 1 < i,j < n, (2.9) 



dA i0 _ dAjo 



0, 1 < i,j < n. (2.10) 



This system may be changed into 

— (—J- 1 ) 
dxi ^ dpj ' dxj \ dp. 



d dJ d dJ 



or 



OXi K Opj ' OXj v opi ' 
after making a transformation 



J- 1 — -)= 0, 1 < i,j < n 



(2.12) 



A 



iO 



9J 7-1 

- — J 1 or A 



,dJ 

OPi 



respectively. The equations (|2.12| ) with n = 2 is a kind of version of the anti-self-dual 
Yang-Mills equations due to Pohlmeyer ]TJ|. Its inverse scattering has been analyzed 



by Beals and Coifman [15] and it may be rewritten in the original Yang's equation "in 



i?-gauge" 16 1 upon choosing 



J 



u 



i w 

v u 2 + vw 



Example 2: Let M — 1, n = 2, ai(A) = —A, 02(A) = A. We obtain the equations 



op 2 dpi 



[A W ,A 21 ] + [A 1U Aw 
[An,A 2 i 



dAn 
dp 2 

dA n | dA 21 _ Q 
9x 2 9xi 



M 21 + M10 M 20 



dpi dx 2 



dxi 



(2.13) 
(2.14) 
(2.15) 



which yields the SDYM equations discussed in Ref. JET], upon making A 2 i — > —A 2 \. 
Here A10, An, v4 2 o, A 2 i are all the Yang-Mills potentials. The corresponding inverse 
scattering problem was introduced by Belavin and Zakharov [|1^] many years ago. 
Recently a scheme of symmetry reduction of the Lax pairs for the SDYM equations 
with respect to an arbitrary subgroup of their conformal group has been described by 
Legare and Popov [T(| and accordingly the compatibility conditions of the reduced Lax 
pairs lead to the SDYM equations reduced under the same symmetry group. 

Example 3: The case of M = 1, Qj(A) = A, 1 < i < n, yields the equations 



[Am, Ajo) — 



dA i0 dAn 



+ 



jo 



dpj dp, 
[A i0) A 3l } + [A a ,A j0 ] - 



0, 1 < i, j < n, 
dA i0 



dA n dA n 
— - H — + 

dpj dpi 



dA 



30 



dxj 



dxi 



[AinAn] + 



dA a dA n 



dxi 



dxj 



0, 1 < i, j < n. 



(2.16) 

0, 1 < i,j < n, (2.17) 
(2.18) 



In comparison with the SDYM equations fl2.13Q ,( [2.14| ), (|2.15|) , these equations may be 
referred to as the generalized SDYM equations. 

Example 4- Let M — 1, a«(A) = A 2 , 1 < i < n. We obtain the equations 



[Aio, Ay 



dA i0 dA j0 
dpj dpi 



0, 1 < i,j < n, 



(2.19) 



dA dA 

[A i0 , A jX \ + [Aa, Ajo] ~ + = 0, 1 < i,j < n, (2.20) 
[A a ,A jX ] + ^ - = 0, 1 < i,j < n, (2.21) 



dAg dAji 
dxj dxi 



0, 1 < i,j < n. (2.22) 



Example 5: The case of M = 1, cij(A) = A + A 2 , 1 < i < n, leads to the equations 

[An, A j0 ] - + = 0, 1 < i, j < n, (2.23) 

dA- dA 

[A i0 , Aji] + [A a , A j0 ] - [A a , A jX \ - + = 0, 1 < i,j < n, (2.24) 

[A iU A n ] + - = 0, 1 < i,j < n, (2.25) 



dAg dAji 

dec j doc <i 



0, 1 < i,j < n. (2.26) 



The systems in Examples 4,5 are new to our knowledge. Obviously each system of 
Examples 3,4,5 is a generalization of Takasaki system of (|2.9| ), (|2.10|) , since all corre- 
sponding reductions of potentials A iX = 0, 1 < i < n, lead to Takasaki system. 

Example 6: Let M > 2, n = 2, a x (A) = -A, a 2 (A) = A M , A X (X) = A 10 + A U X. We 
obtain all equations but the SDYM equations (Example 2) in the SDYM integrable 
hierarchy |]T7| 

[A w ,A 20 ]-^ + ^ = 0, (2.27) 
dp 2 dpi 

i i (Mil cL4 21 <9A 2n 

[A 10 , A 21 ] + L4 n , A 20 - ^1 + ^1 + ^ = 0, 2.28 

dp 2 dpi dx\ 

dA dA 

L4 10 , A 2m ] + [A n , A 2 , m - 1 ] + — ^ + — |^ = 0, 2 < m < M - 1, (2.29) 

[A 10 , A 2M ] + [A n , Vil + ^ + + ^ = 0, (2.30) 

dpi dxi dx 2 

[A U ,A 2M ] + ^ + = 0. (2.31) 

ax 2 dxi 

Note that when M = 2, the equation ( |2.29|) doesn't appear. Recursion operators, 
symmetries and conservation laws associated with this hierarchy have been considered 
by Ablowitz et al ||. In fact, it follows from ( |2.29|) that we have a recursion relation 

A 2m = — (ad^io + ( a ^An + ^-)A 2>m _i = i?y4 2jm _i, 

where R is exactly a recursion operator of the SDYM equations. The above 4-dimensional 
SDYM hierarchy can also be considered as a "universal" integrable hierarchy. Upon 
appropriate reduction and suitable choice of gauge group, it can produce virtually all 
well known hierarchy of soliton equations in 1 + 1 or 1 + 2 dimensions M . 
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3 Darboux transformations 



We recall the spectral problems ( |2.2| ): 

Darboux transformation problem means how we can obtain a new set of 

M 

M\) = J2AiW l<i<n, (3.1) 
1=0 

from a old set of 

M 

V,Ai(\) = ^2A u \ l , 1 <i<n, 

1=0 



so that the spectral problems (|2.2| ) hold covariantly. The purpose of this section is just 
to give an answer to this problem. We shall prove that for certain a to be determined, 

§ = (XI + a)tf, I = diag( 1, 1, • • - , 1 ) (3.2) 

N 

will give a practicable new matrix of eigenfunctions. Of course this matrix a has to 
satisfy some conditions. 

For the time being, let us observe what kind of conditions they should be. Because 
we have 

Li(a)^ = Li(a)[(\I + a)*] = (Li(a)a)* + (A/ + a)L { (a)^ 
= (Li(a)a)*- (A/ + a)i4i(A)* 

r\ M a M 

=Q-T, a ^ k+l ^r) * " ( A/ + «) ( E 4. A 1 < i < ^ 

fc=0 C/Xi i=0 

and 

M 

-i i (A)^ = -(^i,A / )(A/ + «)^, l<^<n, 

a balance of the coefficients of the powers A ; A m , 1 < m < M and A M+1 of Lj(A)\l/ = 
— Ai(X)^f, 1 < % < n, yields that for 1 < i < n, 

da ~ 
aA i0 = -A i0 a, 

OPi 

da 

-Oi,m-i7 Ai m-i - aA im = -Ai >rn -i - A im a, 1 < m < M, 

OXi 

da ~ 

—ttiMT; Am — —AiM- 

OXi 

From these equalities we get the condition for a 

da ~ 

— = aA i0 - A i0 a, 1 < % < n, (3.3) 
dpi 



and at the same time a recursion formula to determine A- TO : 



Am = Am + a iM §^, 1 < i < n, 



A 



i,m—l 



^i,m-l + a i,m-l^: + °A 



zm - A im a, 1 <m < M, 1 < % < n. 



(3.4) 



We observe (3^) and (|3.3|) a little more. Notice that for 1 < r < M, 1 < % < n, we 
may make a further calculation 

A r _i = A r _i + aj , r _i-^ + ai ir + A r (-a) 



A ^ A ( A A \ I \ A 

A,r-i + a i)T „ x - h aA ir +[A ir + a ir - h aA i>r+1 )(-a) + A,r+i( 



M-r 



k=0 

Therefore 



E K 



dxi 

•+fc-l + a i,r+k- 

becomes 



+ «A,r+fc) (-a) fc + A iM (-a) 



-a 



M-r+l 



(3.5) 



4 _ 



k=0 

M-r 



-a 



Ai,r+k—l ^i^+k—lf)^. , 

-!-v/<-< n .A, 1+A (-a) fc , l<r <M, 1< / < /,. 

iM = AiM " / . I / r • 

and further the condition ( |3.3|) reads as 
da 



(3.6) 



Aj 



Am + a iM iz-, 1 <i <n, 



dpi 



aA i0 - A i0 a 



m M da M 

aAo + E Ak(-a) k+1 + aik-^r(-a) h+l + ® A ik(-u) 



k=0 



k=0 



k=l 



M 



M 



da 



M 



\k+l 



a ik o— {-oif 1 " + E aA ik (-a) k , l<i< 

k=0 k=0 aXi k=0 



n. 



(3.7) 



This is a final condition to restrict the matrix a in the construction of new matrices 
of eigenfunctions according to ( ^2|) . The following result can provide us such a kind of 
useful matrices. 

Theorem 3.1 Let h s , 1 < s < N, be N -dimensional column eigenvectors correspond- 
ing to the spectral parameters Ai, A 2 , • ■ ■ , X^, respectively, that is to say that the N- 
dimensional column eigenvectors h s , 1 < s < N, satisfy 



dhs 
dpi 



dh s 

dxi 



- Ai(X s )h s , l<s<N, l<i<n. 



(3.8) 



Assume that the determinant of the matrix H = [hi, h 2 , ■ ■ ■ , hw] is non-zero. Then the 
matrix defined by 

a = -HAH' 1 , A = diag(Xt, A 2 , • ■ ■ , X N ), (3.9) 

satisfies the condition ( \3. Tj) . Therefore we have a Darboux transformation \1/ = (XI + 
a)^ and a new solution An, 1 < % < n, < Z < M, defined by ( \3.dj ). 



7 



Proof: Noting that 



% dy. 
we have for any 1 < i < n 



dH- 1 „ , dH rr , 

- H — H , yi = pi or Xi, 1 < i < n, 



dpi dpi dpi 

OXj OXi OXi 



On the other hand, from (13.81) we obtain 



7T" = E ^fc§^A fe+1 - [A(Ai)/n, • ■ ■ , Ai(X N )h N ) 

OPi k=0 OXi 

M o TT M 

= E a> lk7r A k+1 ~ E AuHA 1 , l<i<n. 

k=(j OX i 1=0 



In this way, we can compute that 



dpi dpi dpi 

M Qtr M 

= ~ E a ik ^-A k+2 H~i + £ AaHA^H- 1 

k=0 OX i 1=0 

M flu M 



+ £ aikHAH-^A^H- 1 - HAH' 1 ]T AuHA'H- 1 

k=0 OXi 1=0 

M fSTT M 

- E a ik ^-H-\-a) k+2 + £ M-<*) W 

k=0 OXl 1=0 

M Qrr M 

+ £ a ik HAH- l ^-H~ l a k + l + a £ ^(-a)' 

fc=0 z=o 

£ ail (-|| Air' + HAH-^H-^(- a f» 

M M 

+ ^A/(-«)' +1 + «EA/(-«y 

M o M M 

fc=0 1=0 1=0 



This shows that the matrix a defined by (|3.9|) satisfies the condition (|3.7|) , indeed. The 
rest of the proof is evident. The proof is completed. I 

We mention that we need a requirement that at least there are two different pa- 
rameters in the set {Ai, A2, • • • , A^} while making Darboux transformation determined 
by Theorem |3.1|. Otherwise we only get an original solution not a new solution since 
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a becomes a unit matrix up to a constant factor. The equation ( |3.8|) for h s is linear 
and hence there is no problem to solve it. So Darboux transformations ( |3.2|) can al- 
ways engender new explicit solutions once one solution is found. Furthermore from the 
solution induced by Darboux transformation, we may make Darboux transformation 
once more and obtain another new solution. This process can be done continually and 
usually it may yield a series of multi-soliton solutions 0, fll3|. 



4 Rational and analytic solutions 



We make Darboux transformation starting from a special initial solution: zero solution 
Au = 0, 1 < i < n, < j < M. Let Ai, A 2 , • ■ ■ , A^v be arbitrary parameter so that at 
least two of them are different. In this way, we obtain a new solution 



M-r+1 
k=0 



da 

h,r+k-l dx 



-a) k , 1< r < M, Ki<n, 



A 



(4.1) 



where the matrix a is defined ( |3.9| ). 

Evidently it is a crucial point to construct an invertible matrix H = [hi, h^, ■ ■ ■ , h^]. 
We shall utilize two special determinants to generate this kind of matrices that we need. 
The one is the Vandermonde determinant formula 



1 

a i 



1 

a-2 



1 

N-l 

N 



Y[(ai - aj), 



(4.2) 



and the other one is a generalized Cauchy determinant formula developed by Constan- 



tinescu IS 



s7V(7V-l)/2 



A 



ii 



A 



lri2 



A nil 



A nirt2 



n 



l -a j ) n 'l\i< l<j <N 2 (b l -b J ) n 



I\l<i<N 1 ,l<j<N 2 ( a i ~~ bj) niri2 



(4.3) 



where niN\ = U2N 2 = N and the matrices A^^, 1 < k < ni, 1 < I < n 2 , are defined 
by 



A, 



'k + l-2 
H '~ 1 fc-1 / , 

' (a^-fti)*^- 1 

When Ni = N 2 = N (this moment rii = n 2 



(ai-6i) fe + i - 1 



(ai-fejv,)^'- 1 



(4.4) 



(ajv 1 -6jv 2 )*+ ! - 1 

1), the above generalized Cauchy 
determinant formula is reduced to the Cauchy determinant formula. 

Let i±, i 2 , ■ ■ ■ , Ziv be N integers and n\Ni = n 2 N 2 = N. For simplicity, we accept 



a(Aj) • p + x = (ai(Aj)pi +x u ---, a n (\i)p n + x n ), 1 < i < N. 



(4.5) 
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We choose the first class of the matrices H as follows 



#1 



ft 
f i 1+ i 



I'l2 

J2 



f. 



2 

12 + 1 



f in 
J N 

J N 



fh+N-1 fi2+N—l e%N 
Jl J2 ' ' ' J N 



(4.6) 



with the functions f\ = /i(a(Aj) ■ p + x), 1 < i < N. This is a little more general 
Vandermonde matrix and hence we have 



detH 1 = ftf?---f N N R(fi-fi)- 



(4.7) 



i>j 



We obtain oji = —HiAiH l , where A x = diag(Ai, A 2 , • • • , A^), provided that detifi ^ 
0. We choose the second class of the matrices H as follows 



Ho 



An(/,0,//) 



Ain 2 (f,g,n) 



&i ni (f,g,iM) ■■■ A nin2 (f,g,fj) 
where Aki(f,g,fi), 1 < k < m, 1 < I < n 2 , are given by 



(4.8) 



Aki(f,g,v) 



'k + l-2" 
, k-l , 



fl 



\l-l)N 2 + l 



_ m )k+l-l 



HN 2 
T N 2 



(^l-fll) 



(/^l-9iV 2 F+^ T 



(;_1)JV 2 + 1 



(MNi-Sl)*^- 1 



/ 



HN 2 
N 2 



(A'iVj —gN 2 ) 



k+l-1 



(4.9) 



with the constants /ij, 1 < i < Ni, and the functions /j = /i(a(Aj) ■ p + x), gi = 
gi(a(\i) ■ p + x), 1 < i < N 2 . The determinant of this matrix H 2 can be computed by 
the generalized Cauchy determinant formula and eventually we obtain 



det# 2 



Y[l<i<N 1 ,l<j<N 2 (f i i — ft')"' 1 ™ 2 

(4.10) 

In this way, we have a2 = — H 2 A 2 H 2 , where A 2 = diag(Ai, • ■ ■ , Atv 2 ; • • • ; Ai, • • •, Aat 2 ) 
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provided that detif 2 7^ 0. 

In what follows, we restrict our consideration within the real field but the case of 
the complex field is completely similar. 

(1) Rational function solutions: 

(1.1) Let fi, ■ ■ ■ , /jv : R n — > R be non-zero distinct rational functions. At this moment, 
we have /j(a(Aj) -p + x) ^ 0, 1 < i < N, and /j(a(Aj) -p + x) ^ fj(a(Xj) -p + x), i ^ j. 
Otherwise we have fi(x) = 0, 1 < i < N, and fi(x) = fj(x), i ^ j, by setting p = 0, 
which contradict to the original hypothesis. It follows from (|4.7[ ) that in this case, detifi 
is a non-zero rational function with the independent variables p, x G E" and hence 
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we can take a special matrix a = a± = —H\K\H 1 , whose elements are all rational 
functions of p, x G R n . Further we can obtain non-zero rational function solutions by 

(1.2) Let fii, 1 < i < Ni, be distinct real numbers, fi : R n —>■ R, 1 < i < N2, be non- 
zero rational functions, and Qi : R n — > R, 1 < i < N2, be non-zero distinct rational 
functions so that Y[i<i<N 1 ,i<j<N 2 (f 1 i — 9j) ^ 0. It follows from fl4.10|) that in this case, 
deti^2 is a non-zero rational function with the independent variables p, x G R n and 
then we can choose a special matrix a = a<i = — H2A2H 2 ~ 1 , whose elements are all 
rational functions of p, x G R n . In this way, we can obtain a class of non-zero rational 
function solutions by ( |4.1|) . 

(2) Analytic function solutions: 

(2.1) Let 71, ■ • • ,7tv be real numbers and hi, 1 < i < N, be any analytic functions to 
satisfy \hi(y)\ < B, 1 < i < N, y G R n . There are a lot of functions of this kind. 
For example, sing(?/), cosq(y), tanhg(y), sechg(y), where q : R n — > R is any analytic 
function. We choose 

fi = fi(a(\i) ■ p + x) = hi(a(\i) -p + x) 1 < i < N. 

This moment, detifi hasn't zero points with respect to (p, x) G R 2n while 

|7«| > B, 1 < i < N, | 7< - 7i | > 2B, 1 < i < j < N. (4.11) 



Therefore under the condition ( 4.11 ) we may take a special matrix a = a\ = —HiA\H 1 1 , 



whose elements are all analytic functions of p, x G R n . This can result in nonzero an- 
alytic function solutions by (|4.1| ). 

(2.2) Let Hi, 1 < i < N ± , be distinct real numbers, fi : R n — > R, 1 < i < N 2 , be 
non-negative or non-positive analytic functions, and hi : R" —> R, 1 < i < N 2 , be 
analytic functions to satisfy \hi(y)\ < B, 1 < i < N 2 . We choose 

Qi = gi(a(\i) ■ p + x) = hi(a(\i) ■ p + x) + 7^ 1 < i < iV 2 , 

where 7 j, 1 < i < N 2 , are all real numbers so that 

|7i - 7j| > 2-B, 1 < « < j < N 2 , | 7i | > B + max |/%|, 1 < i < N 2 . 

l<i<A^i 

In this way, deti^2 hasn't zero points with respect to (p, x) G R 2 ™. Therefore we can 
take a special matrix a = (X2 = —H2A2H2 1 , whose elements are all analytic functions 
of p, x G R n . This may yield a class of nonzero analytic function solutions by means 
of 0. 



5 Discussions and remarks 

From the mathematical point of view, it is very interesting to generate more general 
integrable systems. Such an example may be introduced for the Lax integrable system 
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in Section 2. We display that more general system here 



E [A ik ,Aji]+ J2 ( a i fc 7T^ _ a*fc^M = 0. 

k+l=m k+l=m ^Xj 

-M<k,l<M -M<k,l<M 



dXi 



-2M <m< -M - 1, M + 1 < m < 2M, 



E r , ,in dAi m dAj m 

k+l=m U P] u Pi k+l=m 

-M<k,KM —M<k,KM 



£ ( 



0A, dAu 



%fc-?i Q ifc I — 0, 

— M <m< M. 



where 1 < i, j < n. It is in agreement with the compatibility conditions of the following 
spectral problems with negative powers of the spectral parameter 



( ^ ( E a ' U 

dpi k=—M &3Ci 



M 







M 

-( E A U \> 

l=-M 



1 < i < n, 



where M > 0. A more general Lax set of spectral problems 



9 



M" 



<9 



M" 



E A u \ l p, l<t< 



n. 



l=-M' 



where M', M" > 0, may arise as a reduction of this Lax set with M = max(M', M") 
and vice versa. The Lax integrable system displayed above includes the generalized 
self-dual Yang-Mills flows in Ref. [20|. It seems more reasonable that it is considered 
as a "universal" integrable system, whereas it is too general to lose some concrete 
characteristics. It would be valuable to know whether its more reductions have physical 
interpretations. We may also construct a class of generalized chiral field equations 
similar to one in Ref. 



21 



Essentially the deduction is the same but we should note 
particular properties. 

We remark that there has been a huge class of higher-dimensional nonlinear inte- 
grable equations which can be solved through the nonlocal Riemann-Hilbert method 
proposed by Zakharov and Manakov p2f. This kind of equations are connected with 



the Lax representations defined by higher-order differential operators with respect to 
some independent variables. However, our initial Lax integrable system ( p.5| - |2.8| ) is de- 
rived by the spectral problems involving higher-order powers of the parameter A other 
than higher-order differential operators. Therefore there isn't a direct relation between 
two classes of resulting integrable equations although there exist the same reductions 
of them, for instance, KP equation. On the other hand, our explicit solutions presented 
in Section [| are quite broad because many arbitrary functions are involved in the con- 
struction of solutions and thus it is difficult to discuss their properties. But we should 
be able to study some important properties of sub-classes of the obtained solutions, 
for example, localization property, the existence of the instanton like solutions and 
the soliton interactions etc. as in the work by Zhou f23|. In particular, some careful 



consideration about specific reductions deserves a further investigation. 
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